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Abstract 

The spatial and temporal evolution of small perturbations of the temperature and electromagnetic 
field for the superconducting slab placed in a parallel magnetic field in the regime of thermally 
activated flux creep is studied. On the basis of theoretical analysis of thermal difi'usion and Maxwell's 
equations we find the onset field Bj of the flux-jump instability field and its dependence on the 
external magnetic-field sweep rate Be. 

Key words: thermal and electromagnetic perturbations, critical state, flux creep. 



Introduction 

As wc know that the conventional II type superconductor arc char- 
acterized by a high value of the critical current and upper critical 
field, and is therefore widely used for various technological applica- 
tions. One of the main topic at the technical application of supercon- 
ductors is their stability with respect to flux jumps [1-3]. The flux 
jumps or avalanches are associated with a sudden puncture of mag- 
netic flux into the volume of the sample, and in turn, an increase of 
temperature and the decrease of critical current density. The jump 
phenomena have been observed in conventional hard superconductors 
[1-6], as well as in high-temperature superconductors, recently [7, 8]. 
The critical state stability against flux jumps in hard and composite 
superconductors has been discussed in a number of theoretical and ex- 
perimental papers [1-6]. The general concept of the thermomagnetic 
instabilities in type-II superconductors was developed in literature [4, 
5]. The dynamics of small thermal and electromagnetic perturbations, 
whose development leads to the flux jump, have been investigated the- 
oretically in detail by Mints and Raklimanov [5]. The authors have 
found the stability criterion for the flux jumps in the framework of 
adiabatic and dynamic approximations in the viscous flux flow regime 
of type-II superconductors. Conventionally, thermomagnetic instabil- 
ities were interpreted in terms of thermal runaway triggered by local 
energy dissipation in the sample [5]. According to this theory, any 
local instability causes a small temperature rise, the critical current 
is decreased and magnetic flux moves much easily under the Lorentz 
force. The additional flux movement dissipates more energy further 
increasing temperature. This positive feedback loop may lead to a flux 
jumps in the superconductor sample. 

Theoretical investigations of small thermal and electromagnetic 
perturbations in a various regimes with a various current- voltage char- 
acteristics is one of key problems of electrodynamics of superconduc- 
tors. Moreover, there is a special interest to this problem in the flux 
creep regime of type-II superconductors, since all superconducting de- 
vices used for a large scale applications operate under this regime. In 
our previous work, the dynamics of small thermal and electromagnetic 
perturbations has been studied in the flux flow regime, where voltage 
current-current characteristics of hard superconductor is described by 
linear dependence of j{E) at sufficiently large values of electric field 
E [9]. In the region of weak electric flelds the current- voltage charac- 
teristics j{E) of superconductors is highly nonlinear due to thermally 
activated dissipative flux motion. In the flux creep regime a differential 
conductivity cr strongly depends on the electric fleld E. The nonlinear 
conductivity <t(E) significantly affects the dynamics of thermal and 
electromagnetic processes in superconductors. In particular, it results 
in the dependence of the flux-jump fleld Bj on variation of external 
parameters, in particular on the magnetic fleld sweep rate Bg. A 
theoretical analyze of the flux jumping in the flux creep regime, where 
the current- voltage characteristics of a sample is a nonlinear have been 
carried out recently by Mints [10] and by Mints and Brandt [11]. How- 
ever, a careful study the dynamics of the thermal and electromagnetic 
perturbations in the regime weak electric field with nonlinear current- 



voltage characteristics associated with flux creep is still lacking. 
Objectives 

We report the results of theoretical simulations of the spatial and 
temporal evolution of small perturbations of the temperature and elec- 
tromagnetic fleld for the superconducting slab placed in a parallel mag- 
netic fleld in the regime of thermally activated flux creep. On the basis 
of theoretical analysis of thermal diffusion and Maxwell equations we 
find the onset fleld Bj of the flux-jump instability fleld and its depen- 
dence on the external magnetic- field sweep rate B^. It is assumed that 
the magnetic diffusion is slower than the thermal diffusion. 

1. Formulation of the problem 

Mathematical problem of theoretical study the dynamics of ther- 
mal and electromagnetic perturbations in a superconductor sample 
in the flux creep regime can be formulated on the basis of a system 
nonlinear diffusion-like equations for the thermal and electromagnetic 
fleld perturbations with account nonlinear relationship between the 
fleld and current in superconductor sample. Bean [1] has proposed the 
critical state model which is successfully used to describe magnetic 
properties of type II superconductors. According to this model, the 
distribution of the magnetic flux density B and the transport current 
density j inside a superconductor is given by a solution of the equation 

rotB = j. (1) 
When the penetrated magnetic flux changes with time, an electric 
field E is generated inside the sample according to Faraday's law 

^ dB 

rotE=—. (2) 

at 

The temperature distribution in superconductor is governed by 
the heat conduction diffusion equation 

Kr)^=V[K(T)VT]+J^, (3) 

Here v = y(T) and k = k(T) are the specific heat and thermal con- 
ductivity, respectively. The above equations should be supplemented 
by a current-voltage characteristics of superconductors, which has the 
form 

j = je{T,B)+j{E). 
In order to obtain analytical results of a set Eqs. (l)-(3), we suggest 
that jc is independent on magnetic field induction B and use the Bean 
critical state model jc = jc{Be,T) [1], where Be is the external ap- 
plied magnetic fleld induction. We assume that the dependence of the 
critical current density on temperature of the sample is linear [9] . For 
the sake of simplifying of the calculations, we perform our calculations 
on the assumption of negligibly small heating (T — To <g; Tc — Tq) and 
assume that the temperature profile is a constant within the across 
sample and thermal conductivity k and heat capacity u are indepen- 
dent on the temperature profile; where To and Tc are the equilibrium 
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and critical temperatures of the sample, respectively, [9]. We shall 
study the problem in the framework of a macroscopic approach, in 
which all lengths scales are larger than the flux-line spacing; thus, the 
superconductor is considered as an uniform medium. 

The system of differential equations (l)-(3) should be supple- 
mented by a current- voltage curve j = j{E). In the flux creep regime 
the current-voltage characteristics of type 11 conventional supercon- 
ductors is highly nonlinear due to thermally activated dissipative flux 
motion [12, 13]. For the logarithmic current dependence of the po- 
tential barrier U{j), proposed by [14] the dependence j{E) has the 
form 



-Jc 



(4) 



where the constant parameter n depends on the pinning regimes and 
can vary widely for various types of superconductors. In the case n = 1 
the power-law relation (4) reduces to Ohm's law, describing the normal 
or flux- flow regime [15]. For infinitely large n, the equation describes 
the Bean critical state model j = jc [!]■ When 1 < n < oo, the equa- 
tion (4) describes nonlinear flux creep. In this case the differential 
conductivity a is determined by the following expression 



3. Dispersion relation 

Let us derive the dispersion equation to determine the eigenvalue 
problem. As we know [5] that a nature of the flux jumps depends on 
the competition between diffusive and dissipative processes through 
the dimensionless parameter 
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where Dt = nfv is the thermal diffusivity and Dm = c? /i-Kaf the mag- 
netic diffusivity coefficients, respectively. Therefore the flux instabil- 
ity criterion is determined mainly by the relation of the magnetic Dm 
and thermal Dt diffusion coefficients. As we have mentioned above, 
the differential conductivity <t(E), which determines the dynamics of 
the instability is high in the flux creep regime and the parameter t 
is high enough, also. It is clear that this picture for the flux jumps 
corresponds to the limiting case r 2> 1- Consequently, it can be as- 
sumed that the initial rapid heating stage of a flux jump takes place on 
the background of a "frozen-in" magnetic flux. Therefore, under this 
dynamic approximation, we obtain from (9) the relation between elec- 
tric field t{x, t) and temperature Q{x, t) perturbations in the following 
form 



(5) 
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According to relation (5) the differential conductivity decreases with 
the increasing of the background electric field E. and strongly depends 
on the external magnetic field sweep rate Ei, ~ B^x. Therefore the 
stability criterion also strongly depends on the differential conductiv- 
ity a. For the typical values of ji = 10**A/cm'^, Ei, = lO^^V/cm 
we obtain a = lO^^l/ficm. It follows from this estimation [11] that 
the differential conductivity a which determines the dynamics of ther- 
momagnetie instability for is high enough. We assume, for simplicity, 
that the value of n temperature and magnetic-field independent. 

2. Basic equations 

Let us formulate a differential equations governing the dynamics of 
small temperature and electromagnetic field perturbation in a super- 
conductor sample. We study the evolution of the thermal and electro- 
magnetic penetration process in a simple geometry - superconducting 
semi-infinitive sample a; > 0. We assume that the external magnetic 
field induction Eg is parallel to the z-axis and the magnetic field sweep 
rate Eg is constant. When the magnetic field with the flux density Eg 
is applied in the direction of the z-axis, the transport current Sj{x,t) 
and the electric field SE{x, t) are induced inside the slab along the 
y-axis. For this geometry the spatial and temporal evolution of small 
thermal ST{x,t) and electromagnetic field 5E{x,t) perturbations 

5T = 0(a;)exp[7t], (6) 

= e(a;)exp[7t]. (7) 

is described by Maxwell equations coupled to the thermal diffusion 
equation 
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where 7 is the eigenvalue of the problem to be determined. It is clear 
that the rate 7 characterizes the time development of the instability. 
In the case when Re7 > 0, small thermal and electromagnetic pertur- 
bations increase and the stability margin corresponds to the case when 
7=0. It should be noted that the nonlinear diffusion-type equations 
(8) and (9), totally determine the problem of the space-time distribu- 
tion of the temperature and electromagnetic field profiles in the fiux 
creep regime with a nonlinear current-voltage characteristics in the 
semi-infinite sample. 
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We notice that the last relation between e(x, t) and temperature &{x, t) 
has been derived in the assumption that the decrease of the critical 
current density jc resulting from a temperature perturbation &{x, t) 
compensates with increase of the resistive current density jr resulting 
from an electric field perturbation e(x,t), so the total current density 
remains constant [5]. Upon substituting the expression (10) into the 
equation (8) and excluding the variable e{x,t) one can get the differ- 
ential equation for the distribution of thermal perturbation, which can 
be conveniently presented in the following dimensionless form 



dp^ 



- p0 = 0. 



(11) 



Here we introduced the following dimensionless variables 
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Here L 



cBe 



is the magnetic field penetration depth. Thus, the 



condition of existence of a non-trivial solutions of Eq. (11) allows to 
define the spectrum of eigenvalues of 7 and the instability threshold, 
accordingly. The equation (11) has an exact solution in terms of Airy 
functions given as the following form 

0(p) = ci(s)Ai(p) + C2{s)Bi{p). (12) 
where Ai(p) and Bi[p) are the Airy functions. Here constants of in- 
tegration ci and C2 are determined from the thermal boundary condi- 
tions. Substituting the last solution (12) into the thermal boundary 
conditions (see, for example [9]) we find that C2 = and Q{p) = 
ci{s)Ai{p). Applying the second boundary condition 0(1) = we get 
an equation to determine the eigenvalues of the problem 

where an are the zeros of the Bessel function and growing with in- 
creasing n; For example, for n=l the stability criterion is presented 
as 

01 = r^/^-y. (13) 
Using the value for the magnetic field penetration depth, we can easily 
obtain from (13) an expression for the threshold magnetic field Bj at 
which the branching instability occurs 
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Let us now estimate the threshold field for a typical values of param- 
eters jc ~ lO^A/m^, Tc - To ~ 10 K, K ~ IQ-^W/Km, n=10. The 
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background electric field Ej^ ~ B^L, induced by the magnetic-field 
variation Be ~ lO^^ ^ 10'^ T/s is of the order of = 10~^ ^ lO'^ 
V/m for the value of L = 0.01 m. We can easily estimate that the 
threshold field has the value X 1 3T. 

4. Discussion 

Experimentally, the background electric field is created by the 
sweeping rate of applied magnetic field Be- As can be seen from the 
relation (14) the threshold field Bj is decreased with the increasing of 
background electric field. It is noticeable that the dependence of the 
flux-jump field Bj on the sweeping rate Be of the applied magnetic field 
have been verified by a numerous experiments [8, 16-18]. An inten- 
sive numerical analysis on the sweep rate dependence of the threshold 
field has been performed recently in [19]. Recent magnetization mea- 
surements [8, 16] have shown that the value of the threshold field Bj 
decreases as the sweep rate Be increases. A theoretical investigations 
on the dependence of threshold field on the varying external magnetic 
field has been performed in detail recently by Mints [10]. Within the 
framework of the flux jump instability theory [4, 5] a rapid variation 
of the applied magnetic field acts as the instability-driving perturba- 
tion, and that threshold field Bj should decrease with increasing the 
sweeping rate Be [8]. The numerical studies [19] have demonstrated 
that the flux jumps takes place when the sweep rate Be increases up to 
a certain value, where the number of jumps increases with the sweep 
rate Be- As the sweep rate further increases, these simulation results 
show that the flux-jump fleld decreases and approaches a saturation 
value, which is fairly close to the experimental value of about 1.2 T/s 
[8]. However, as has been mentioned in [16], an experimental investi- 
gations on the dependence of the threshold fleld for the flux jump fleld 
Bj on the external magnetic field sweep rate Be is very little. Exper- 
imentally, can be observed a complex behavior of dependence of the 
threshold field Bj on the sweep rate Be ■ The results of experiments of 
Ref. [16] demonstrated that Bj is independent of the sweep rate in a 
defined range of temperatures. Thus, the near independence of Bj on 
the sweeping rate remains to be explained. In some conventional su- 
perconductors both the independence of Bj on the sweeping rate [3, 4] 
and its growth at a high sweeping rate [20] were detected. It has been 
suggested [4] that a nonuniform heating may be responsible for such 
an effect. However, theoretical understanding of the thermomagnetic 
instabilities at such conditions is still lacking. Note, however, that 
some details of the local fleld behavior depend indeed on the sweeping 
rate, as, for example, the number and amplitude of the jumps. Gerber 
at. al. [18] have demonstrated that at low values of the sweep rates 
the number of flux jumps decreases as sweep rate increases. At still 
high sweep rates the amplitude of flux jumps becomes independent of 
the sweep rate and saturates to the limit with further increasing sweep 
rate. 
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Fig.l The sweep rate dependence of the fleld of the flux jump field. 



In Fig.l we have demonstrated the dependence of the threshold 
fleld Bj on the external magnetic field sweep rate. As can be seen, 
the value of Bj decreases as the sweep rate increases. As the sweep 
rate increases the value of Bj decreases and it tends to saturate at 
high sweep rates [8] . Magnetic fleld dependence of the critical current 
density only slows down the decrease of the field Bj with increasing 
external magnetic fleld sweep rate [10, 11]. We note that for the case 
Kim- Anderson model [12], the absolute value of the exponent in the 

■ —1/3 

power formula decreases from 1/2 to 1/3, so Bj ~ Be [10]. 

Let us qualitatively estimate the temperature dependence of the 
flrst flux jump fleld Bj{T). In order to compare the last formula with 
experimentally determined flrst flux jump fleld, one must know the 
temperature dependencies of the critical current density and the spe- 
cific heat of the sample. As has been mentioned in literature [5, 16] a 
quantitative estimation of ic(T) from the available experimental and 
theoretical models is difficult because of the uncertainty in the values 
of the critical current jc at a given field and temperature. To deter- 
mine jc{T), different approaches have been taken. Empirically, the 
critical current density jc(T) can be presented in the form 

jc = jo [1 - t"]™ . 

where 1 < m < 2, 4 = T/Tc- The different exponents n=l and 2 
refer to the most common cases discussed in the literature, where the 
critical current exhibits a linear and a quadratic dependence on T/Te- 
There is experimental evidence [21], which indicates that the temper- 
ature dependence of the critical current density approximately linear 
at low temperatures. However, at higher temperatures, where flux 
creep effects are dominant, the temperature dependence of the critical 
current density can be presented as 

ic(T) = >(o)(l-^2)^ 

commonly accepted in literature [6]. Assuming that the thermal con- 
ductivity is a linear function of temperature, we can easily obtain an 
expression for the temperature dependence of the threshold fleld Bj (T) 
(Fig 2.) 

BjiT)^^t{l-t'2f. 
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Fig. 2. The temperature dependence of the flux jump field 
Conclusion 

We have performed a theoretical study of dynamics of small ther- 
mal and electromagnetic perturbations in type-II superconductors in 
the flux creep regime with a nonlinear current-voltage characteristics. 
For this purpose, the space-time evolution of temperature and electric 
fleld was calculated using the heat diffusion equation, coupled with 
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Maxwell's equations and material law, assuming that the ajiplied mag- 
netic field is directed parallel to the surface of the sample. We found 
the threshold field Bj for the occurrence of thermomagnctic instability 
assuming that the heat flux diffusion is considerable faster than the 
magnetic flux diffusion in superconductor. The obtained stability cri- 
terion for the thermomagnetic flux jumps demonstrates the extremely 
high sensitivity of the threshold field Bj on the values of the criti- 
cal current density jc, thermal conductivity «, and external magnetic 
field sweep rate Be- Thus, the stability condition (14) for the ther- 
momagnetic flux jumps directly reflects the magnetic sweep rate Be 
dependence on the threshold fleld Bj. It follows from the criterion 
that the value of the threshold field Bj is inversely proportional to the 
square root of the magnetic-field sweeping rate Be. Therefore, with 
the increase of sweeping rate Be the threshold field Bj decreases. Fi- 
nally, we have discussed the temperature dependence of the threshold 
field, briefly. 
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